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On the possibility of using seismic probes to study the core composition 
in pulsating white dwarfs 
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White dwarfs correspond to the final stages of stellar evolution of solar-type stars. In these objects, production of energy 
by nuclear burning has ended which means that a white dwarf simply cools down over the course of the next billion years. 
It is now known that white dwarfs spend some of their cooling history in an instability strip. The pulsating white dwarfs 
with an hydrogen atmosphere (called DAV or ZZ Ceti stars) show non-radial oscillation modes with periods in the range 
100 - 1200s. In this work we try to illustrate how the oscillation p-mode frequencies of idealized white dwarf models 
change as the result of a different chemical composition in the core, with the ultimate goal of determining the chemical 
stratification from seismic observations. The presence of acoustic glitches in the internal structure results in a periodic 
signal in the frequencies. We find that this signal depends on the chemical stratification/composition of the core in a form 
that can be analytically modelled. 
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1 Introduction 

At the end of the main sequence phase, a small mass star 
has exhausted the nuclear fuel in its core. It looses a signif- 
icant part of the mass of its outer layers while in the AGB 
phase. Depending on when this mass loss ends (by dissipa- 
tion of the rarefied expanded envelope) the carbon-oxygen 
core of the star will be surrounded with either a thin coating 
of hydrogen-rich material or a layer of nearly pure helium. 
We call the remaining star a white dwarf. 

The hydrogen atmosphere (DA) white dwarfs comprise 
about 80% of all spectroscopically known white dwarfs and 
are important, e.g., in u nderstanding the mass loss process 
in the AGB phase (e.g. lMarig o 2012) and determining the 
age of the local galactic disk dWoodlll992[). Since the first 
detection of a pulsating white dwarf dLandolt 1968 ) the 
number of known pulsators has grown to about two hundred 
and it is now believed that these stars go through instability 
strips during the cooling process. The instabilities manifest 
in terms of multi-periodic luminosity variations which can 
be explained by non-radial g-mode pulsations. This means 
that we can use the tools of asteroseismology to probe the 
internal structure of white dwarfs and determine some of 



their fundamental properties (e.g. Bradley & Winget 1994 



Cdrsico et al. 2012). For a review on the recent state of 



seismo lo gy studies on white dwarfs see IWinget & Keplei 
d2008l) or lAlthaus et all d2010h . 



Despite never having been detected observationally 



dSilvotti et alJ 1201 II) . pressure-driven modes (p modes) 
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could also be excited in white dwar fs. These are exp ected to 
show periods between 0.1 and 10 s ( Saio et alll983 ). These 
modes probe the interior of the star differently as they are 
more sensitive to the structure of the core (g modes are en- 
velope modes). Our goal in this work is to study the effect 
that a different chemical composition in the core has on the 
oscillation frequencies of p modes. 

Very simple models of the structure of white dwarfs are 
computed using the appropriate equations (this is discussed 
in $2}. The sharp transition caused by a chemical compo- 
sition discontinuity gives rise to a characteristic signature 
in the oscillation frequencies. The nature of this signature, 
in the case of p modes, is discussed in |j3]and in fj4] results 
from the idealized models show the core composition de- 
pendence of the signal. The paper ends with a discussion of 
the results. 

2 Equilibrium Models 
2.1 The Equation of State 

The pressure associated with matter at T=0 is what keeps 
the star from collapsing. The equation of state of this matter 
can be seen as that of an ideal, cold gas of fermions. Follow- 
ing a usual convention we define the dimensionless Fermi 
momentum, usually called relativity parameter, by x=-^r- 



so tha t we can write the pressure of the gas as (Weis s et al 
2004) 



P : 



87rm|c 5 
3^ 3 



,.'4 



(x' 2 + 



-dx' = Af(x) (1) 
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and the mass density 



5.0 



3h 3 N, 



(J, e X 



B p e x i 



(2) 



where A and B are constants, m e is the mass of the 
electron, h the Planck constant, c the speed of light and Na 
the Avogadro number. The symbol p e represents the inverse 
of the number of electrons per unit atomic mass, sometimes 
refered to as the mean electron weight. 

The essential role of the equation of state in this case is 
to specify the mechanical structure of the models, to which 
pulsation modes are very sensitive. Because the pressure in 
the interior of white dwarfs is completely dominated by the 
contribution of the cold and highly degenerate electron gas, 
we can expect this equation of state to be adequate in con- 
structing models of pulsating white dwarfs. The goal here 
is not to develop a realistic model of the structure of white 
dwarfs but one that includes only the necessary physics to 
yield reliable frequency calculations. 

2.2 The Chemical Composition 

We have built several white dwarf models by integrating the 
hydrostatic equilibrium equation together with the equation 
of state defined by (Q~|) and (f2]). The value of the central den- 
sity is provided as input and iterated until the model has 
mass M, 

We consider an inner core of a specified chemical com- 
position and mass surrounded by an envelope of pure hy- 
drogen, thus neglecting the He layer that should be present 
in DA-type stars. The composition discontinuity in the core- 
envelope transition is handled by changing the value of the 
mean electron weight fi e . 

The cores of typical white dwarfs are expected to 
be composed of the products of He burning: carbon 
and oxygen. There is however some observational and 
theoretical basis to expect up to iron-r ich cores (see 
lHamada & Salpeterlll96ll:|Panei et alJboOOl) . We define the 
following name convention for our models: model Mj has 
a core composition of either j = C (Carbon), O (Oxygen), I 
(Iron), M (Magnesium). We have also considered homoge- 
neous Carbon models without the hydrogen layer, for com- 
parison. 

3 Effect on the frequencies 

The discontinuity in the chemical composition at the core- 
envelope interface is expected to have an impact on the pul- 
sation frequencies of the models, mainly due to its influence 
on the Brunt-Vaissala frequency N, which is defined as 



N 2 = g 



1 d\nP d\np 
Ti dr dr 



(3) 



where g=Gm/r 2 is the gravitational acceleration, G the 
gravitational constant and Ti=(d\nP/d\np) s is the first 
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Fig. 1 The square of the Brunt-Vaissala frequency (solid 
lines) and the square of the sound speed (dashed lines) as 
a function of the natural logarithm of the ratio of radius to 
pressure. This coordinate is used to provide appropriate res- 
olution at both the center and the surface of the model. The 
two colors represent models with different masses: 0.8 M Q 
in red and 0.9 M Q in blue. The two models have slightly 
different radii. 

adiabatic exponent (derivative at constant entropy s). This 
is the natural characteristic frequency that dominates the 
calculation of pulsation periods in white dwarfs. Any sharp 
feature (in our case a discontinuity) on the run of as a 
function of depth (see Figure [TJ will have an effect on the 
frequencies. In this section we briefly describe a method 
used in solar-type stars that allows us to predict what this 
effect will be, in the case of p modes. 

Following dChristensen-Dalsgaard et al .1 19951) . let V(t) 
be the acoustic potential which can be written as 



V = N + 

c d 
~2~dr 
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-4nGp, 



where c is the sound speed and r is the acoustic depth 
defined by 



r(r) 



dr' 

c 



(5) 



The acoustic size of the star is r s = t(0), while R is 
the radius of the star. A discontinuity in the sound speed 
(Figure [TJ due to the chemical stratification gives rise to a 
discontinuity in the potential V(r). 

The eigenfunctions Y of the modes of oscillation sat- 
isfy, asymptotically and f or low degree t, the equation (e.g. 
Vorontsov. Zharkovl 19891) 



d 2 Y 
dr 2 



+ [io 2 - V 2 (t)] Y = 0, 



(6) 
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Fig. 2 Acoustic potential as a function of the acoustic 
depth r (note that r = is the surface). The four models 
shown have 0.9M© and different core compositions. The 
black line corresponds to the homogeneous carbon model, 
without the hydrogen layer. The step in the potential has a 
different amplitude and location for each model. 



where uj is the frequency of the mode. For the case of 
a potential with a step function of amplitude 5V, we can 
solve equation §6^ an d show that th ere will be a periodi c 
signal arising in the frequencies (see lMonteiro et al.lll994l) . 
The change in the the frequencies, 5uj, relative to the unper- 
turbed frequency uiq, will be 



SV 2 

5uj ~ 7 sin (2ujQT d ) = A 



sin(2w Trf) , (7) 



with Td the acoustic depth at which the discontinuity oc- 
curs. We use a reference frequency w r /2ir = 10 Hz (repre- 
sentative value within the frequency range for p modes). 

Note that both the amplitude (which is a decreasing 
function of frequency) and the periodicity of the signal are 
completely determined by the amplitude and location of the 
step in the potential. In Figure ©, V(t) is plotted for four 
different models, each with a different core composition. 
For the homogeneous model there is no step function and 
thus we expect no periodic signal. 

4 Results 

For models with 0.8 M© and 0.9 M©, the results for the am- 
plitude Aq and the acoustic depth r s of the discontinuity on 
the potential are shown in Figure (O. Each different com- 
position results in a different point in this graph, meaning 
that the signal in the frequencies will depend on the core 
composition of the model. 

Using the POSC code dMonteiroll2008l) . linear adiabatic 
frequencies were calculated for each model, for modes with 
I = 0, 1 and 2. The large separations, defined as the differ- 
ence between the frequencies of consecutive modes of the 
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Fig. 3 Results for the amplitudes and acoustic locations of 
the discontinuity for models with 0.8M© (red) and 0.9M© 
(blue). Each point represents a different core composition, 
as indicated by the labels. 




Fig. 4 Large separation as a function of frequency for 
three representative models of 0.8M© with an homogeneous 
carbon composition (solid line), a carbon core (dotted line) 
and an iron core (dashed line). Filled circles represent I — 
modes, open circles £ = 1 and stars 1 = 2. 



same degree £, were also computed. The predicted periodic 
signal appears on this quantity, as is visible on Figure ©, 
for three representative models of 0.8M©. The trend in the 
amplitude of the signal with the frequency of the mode is 
noticeable in the Mr model. The periodicity is also differ- 
ent for the two models with a discontinuity. 

5 Discussion 

Very simple equilibrium models of white dwarf stars were 
built in order to calculate oscillation frequencies. The hy- 
drostatic equilibrium equation was integrated, considering 
the equation of state of a degenerate, ideal gas of electrons. 
It is clear that this part of our work requires a more ade- 
quate approach and that a study with more realistic models 
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(full evolutionary models that take into account the cooling 
process and a more realistic equation of state) is required. 

The equation of state is not the major problem though: 
the treatment of the chemical stratification done here is sim- 
plistic. However, our models contain the relevant aspects of 
the physics of cool white dwarfs that may be probed with 
seismic data and are, thus, capable of providing useful re- 
sults. 

Calculating the size and location of the step function in 
the acoustic potential V{t) allowed us to predict the effect 
that the core-envelope discontinuity has on the p-mode fre- 
quencies: a periodic signal. This signal is clearly seen in the 
model frequencies (Figure |4|i and depends substantially on 
the chemical composition of the core. The absence of any 
periodic variation in the homogeneous model also confirms 
our predictions. The periodic signal in the large frequency 
separation displays a variation of amplitude with frequency 
and a periodicity associated with the acoustic location of the 
discontinuity, as predicted. These results open the window 
for further detailed analysis of both more accurate models 
and observational seismic data of white dwarfs. 

The asymptotic analysis considered in !j3] is valid for 
p-mode frequencies and is commonly used in the context 
of heliose ismology and asteroseismology of solar-type stars 
(see e.g jMonteiro et alJl200(A iMazumdar et al.ll2012h . 

We emphasize again that, until now, only g-mode oscil- 
lation frequencies have been detected in white dwarfs. The 
very high frequencies expected for p modes make these very 
hard to detect with present day instruments, even if they are 
indeed excited. This does not invalidate our analysis as the 
sharp features also gi ve rise to a similar ef fect in the period 
spacings of g modes dMetcalfe et al ] l200ll) . 
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